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We study transport properties of a non-degenerate two-dimensional system of non-interacting 
electrons in the presence of a quantizing magnetic field and a short-range disorder potential. We 
show that the low-frequency magnetoconductivity displays a strongly asymmetric peak at a nonzero 
frequency. The shape of the peak is restored from the calculated 14 spectral moments, the asymp- 
totic form of its high-frequency tail, and the scaling behavior of the conductivity for uj ^ Q. We also 
calculate 10 spectral moments of the cyclotron resonance absorption peak and restore the correspond- 
ing (non-singular) frequency dependence using the continuous fraction expansion. Both expansions 
converge rapidly with increasing number of included moments, and give numerically accurate re- 
sults throughout the region of interest. We discuss the possibility of experimental observation of the 
predicted effects for electrons on helium. 



I. INTRODUCTION 



Single-electron dynamics in the lowest Landau level 
(LLL) broadened by a delta-correlated scalar disorder 
potential provides the simplest framework for the analy- 
sis of the integer jnuantum Hall effect (IQHE). Wegner's 
exact calculationld of the density of states raised hopes 
that the corresponding model may be exactly solvable, 
and much effort was put into understanding transport in 
this model. However, in contrast to the density of states, 
the conductivity is expressed in terms of a two-particle 
Green's function and depends not only on the energies of 
single-particle states, but also on their wave functions. 

A noteworthy feature of the single-electron wave func- 
tions in a strong magnetic field is delocalization near the 
centers of the disorder-broadened Landau bands. For 
small energies E counted off from the band center (in 
the neglect of band mixing) ,r\ke localization length ^ di- 
verges as a universal powedj'cl of \E\. These are tran- 
sitions between large-radius states that form the low- 
frequency conductivity of the system. As a result, for 
Fermi energies close to a band center, the model displays 
a universal critical behavior at sufficiently small temper- 
atures and frequencies. 

The width of the critical region depends on the proper- 
ties of the disorder potential and the Landau level num- 
ber. For the lowest Landau level, it is of the order of 
the bandwidth ?i7, the only dimensional parameter of 
the Hamiltonian projected on the LLL (we assume that 
the cyclotron frequency lUc ^ 7)- Outside the scaHng 
region, the spatial extent of the eigenstates is small, of 
the order of the magnetic length I = {h/muscy'/'^ , and 
the universality is lost. Therefore the overall frequency 
dependence of the conductivity is determined by the dis- 
order mechanism and may allow to discriminate between 



different mechanisms. For this reason, it is interesting 
to obtain the frequency dependence of the conductivity, 
including its universal and non-universal parts, at least 
for some basic models of disorder. It is also of interest 
to find highly accurate numerical results, as they may be 
used to test various approximate analytical approaches. 

In the present paper we consider the frequency- 
dependent conductivity for a short-range disorder poten- 
tial. Short-range disorder is usually a good model of the 
random potential experienced by electrons trapped on 
the surface of liquid helium. This potential is due to qua- 
sistatic ripplons, slow helium ijapor atoms, or, for thin he- 
lium films, substrate disordero. The electron correlations 
in this system are usually strong. However, it was be- 
lieved that, for strong enough magnetic fields, the static 
magnetoconductivity can be described in the single- 
electron approximation (for a review see Refs. |-0). The 
corresponding results obtained within the self-consistent 
Born approximations (SCBA) appeared, to be in reason- 
able agreement with the experimentETO. 

The SCBA ignores the interference effects which lead 
to electron localization and to the scaling behavior near 
the band center. Therefore there is an apparent contra- 
diction between the interpretation of the experimental 
data on electrons on helium and the phenomenology of 
the quantum Hall effect. From this point of view, it is 
important to develop a single-electron theory of magne- 
totransport that will take localization effects into account 
and will thus extend the ideas of the IQHE theory to a 
new parameter range, which is of interest for electrons on 
helium in particular. The results can serve as a basis for 
the full many-electron-t|heory of magnetoconductivity at 
strong magnetic fieldalj. 

With this in mind, we will analyze magnetotransport 
in the regime 
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h-i < ksT, < 1, (1) 

where the disorder-induced broadening of the lowest Lan- 
dau level is small compared to temperature, and simul- 
taneously the filling fi'action is small, i.e., the electron 
system is nondegenerate. Much of the experimental data 
on electrons on helium refers to the range (Q) , as the sys- 
tem is very clean, and the experiments are often done at 
low electron densities n ~ 10* cm~^. This range is also 
of interest for low-density electron systems in semicon- 
ductors. 

In the range all states within a broadened Landau 
level are nearly equally populated, and one no longer 
needs to take into account the Boltzmann factor while 
computing temperature and disorder averages using the 
Hamiltonian projected on that level. As we show below, 
this simplification allows us to make accurate calcula- 
tions of the frequency-dependent conductivity using the 
method of moments, which, was previously suggested for 
this problem by one of usllj. 




FIG. 1. Reduced microwave conductivity of a 

non-interacting 2DES at the lowest Landau level in a 
short-range disorder potential for lj,7 ^ ksT/h. For small 
frequencies, a; <C 7, the conductivity is singular, a^x ~ i-o'^. 
It is determined by spatially large, nearly delocalized states. 
For oj ^ 7, the conductivity is determined by large optimal 
fluctuations of the disorder potential as illustrated in the in- 
set. The optimal potential V'opt(r) is such that hu; equals to 
the energy difference Et — Eh between the top and bottom 
bound states \t), \b), and at the same time these states are 
maximally overlapping. 

The outline of the paper is as follows. In Sec. II we cal- 
culate the frequency-dependent conductivity axxi^) in 
the temperature range h'^ <C ksT, exp{huJc/kBT) ^ 1 
for uj LUc [these results were previously announced 
in Ref. |l^. We find the asymptotics of the conduc- 
tivity at both small (a> — > 0) and comparatively large 
(7 ^ cj <C LOc) frequencies, and show that axxi^^) has 
a peak at a nonzero frequency a; ~ 7. Using an effi- 
cient diagram classification scheme, we compute exactly 



the first 14 spectral moments of this peak. These mo- 
ments contain information about the short-time (~ 7^^) 
dynamics of the system. Combined with the low and 
high-frequency asymptotics, they allow us to accurately 
restore the entire function axxi^) [see Fig. In Sec. Ill 
we investigate the cyclotron resonance, i.e., axxii-^) for 
w ~ Wc- We calculate the first 10 frequency moments of 
the cyclotron resonance absorption peak and use them to 
accurately restore its shape [see Fig. |j as a function of 
frequency detuning Auj = oj — ujc- To do the restoration, 
we also calculate the asymptotic form of the tails of the 
cyclotron resonance using the method of optimal fiuctua- 
tion [our result differs from that obtained earlier by loffe 
and LarkinEJ]. In Sec. IV we discuss the ways to observe 
the predicted here behavior in experiment. Technical de- 
tails are given in the Appendices. 
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FIG. 2. Reduced conductivity at the cyclotron absorption 
peak [see Eq. (^)] of a non-interacting 2DES in a short-range 
disorder potential for \uj — uic\, 7 <C ksT /% <^ uj^. The curve 
is analytic in the center of the peak. The absorption at the 
tails, — LJc| 3> 7, is determined by large optimal fluctuations 
of the disorder potential as illustrated in the inset. 



II. LOW-FREQUENCY SINGLE-ELECTRON 
CONDUCTIVITY 

In this section, we calculate the conductivity of a non- 
degenerate non-interacting two-dimensional electron sys- 
tem (2DES) for low frequencies, uj <^ ksT/h. We con- 
sider the case of a delta-correlated disorder potential 
and assume that the broadening of the Landau level 
7<u;e, kBT/h jcL Eq. @)]. 

In the range d3) all states within the LLL are equally 
occupied, and the Kubo formula for the dissipative con- 
ductivity can be written as a simple trace without the 
Boltzmann factor, 

Ox. {^) = -^—^ 3f^e dt e'"' {jx it) jx (0) ) (2) 
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Here jx = epx/m is the one-electron current operator, 
/3 = Ti/ksT, the angular brackets (•) denote statistical 
averaging over the states followed by an averaging over 
quenched disorder, while the horizontal line denotes only 
the disorder averaging. The trace Tro in Eq. (|) IS per- 
formed over all single-particle states of the lowest Landau 
level- the energies are measured with respect to its center. 
Eq. (0) is written for the case of strongly quantizing mag- 
netic fields, exp(/3a'c) ^ 1, so that only the lowest Lan- 
dau level is occupied. However, the calculation is readily 
generalized to the case of arbitrary l3ujc by replacing Trg 
by the sum of traces over the states of each Landau level 
n weighted with exp(— n/3tt'c)[l — exp(— /?cjc)]- 

Calculations within a single Landau level are conve- 
niently done using the formalism of the guiding center 
coordinates R = {X,Y). The electron dynamics in the 
random potential V(r) is mapped onto that of a ID quan- 
tum particle with the generalized momentum and coor- 
dinate X and Y, and with the Hamiltonian 

H ^h-fS^ V„ exp(iqR), [X, Y] = -if. (4) 

Respectively, the guiding center velocity is determined 
by the potential gradient, 



where D = P'j is the characteristic diffusion coefficient 
and, as discussed in Appendix 



die*"*Tro{R(t) -11(0)} 



(10) 



is the reduced conductivity. It depends on the ratio uj/"f 
of the only two quantities with the dimension of fre- 
quency that remain after projection on one Landau level. 

The expression (|o|) can be rewritten with the help of 
Eq. d) as 

/OO 
die-*^ ,(qq') 

xTto{%\' exp [zqR(0] exp [^q'R(0)]}. (11) 

This form is particularly convenient for calculating the 
frequency mom ents of the reduced conductivity, see be- 
low in Sec. HC . 

Yet another representation of the reduced low- 
frequency conductivity can be obtained if we describe 
time evolution of the electron operators in Eq. ( pi] ) using 
the set \n) of the eigenstates of the full electron Hamil- 
tonian for the lowest Landau level and perform the time 
integration. 
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Y,KEn-Ern-nu;)\{n\VV\m)\ , (12) 



where fi,!/ — x, y, and e^^ is the unit antisymmetric ten- 
sor, f^xy ^yx — 1- 

The dimcnsionlcss coefficients 



Vq = (VMexp(-;V/4) 



(6) 



are proportional to the Fourier components of the disor- 
der potential. 



(7) 



where S is the area of the system. For higher Landau lev- 
els the coefficients have to be modified as explained 
in Appendix ^ [see Eq. (A3)]. We will assume that V^(r) 
is zero-mean Gaussian and delta-correlated, 



V{r)V{r') = v^S{r-r') 



(8) 



in-which case the SCBA width of the lowest Landau band 
isEl?i7= {2/ny/^v/l. 

In the simplified Kubo formula (|^) the temperature 
dependence is factorized, and we can rewrite the low- 
frequency conductivity in the form of the generalized Ein- 
stein relation 



ne'^D 1 



(9) 



where £"„ are the energies of the LLL states \n) in the po- 
tential V{r) (again, generalization to the case of several 
occupied Landau levels is straightforward). 

We emphasize that, in the chosen parameter range, 
the Landau-level projection resulted in expressions that 
do not contain the usual disorder-dependent denomina- 
tor, and the quenched disorder averaging can be done 
directly, without invoking supersymmetry or the replica 
trick. 



A. Tail of the low-frequency conductivity 

We begin with calculating the asymptotic form of the 
reduced conductivity (T(a;) for w ^ 7 from Eq. (^2|). In 
the neglect of inter-band mixing, the energies En are 
symmetrically distributed around the Landau band cen- 
ter {E ~ 0). The tails of the density of states p{E) are 
known to be Gaussian, p{E) cx exp(— 4i5^/?i^7^). They 
are determined by the probability of the optimal (least 
improbable) potential fluctuation Ve (r) in which tlie. low.- 
est or highest bound state has energy E (|^S> 7)liJ'BEj. 

If we ignore the matrix element in Eq. (|lj) altogether 
(as we show below, this only affects the prefactor), the 
tail of the conductivity will be proportional to the proba- 
bility to find two states En, Em such that En — Em = huj- 
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The major contribution comes from states at the op- 
posite ends of the energy band with energies close to 
En = -E,n^'huj/'2., giving 
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[p{riuj/2)Y oc exp(-2w77^) 



(13) 



To check this approxi«i^tion, we will apply the method 
of optimal fluctuationtSEJ. The averaging over disorder 
in Eq. ( p2| ) will be done using the path integral represen- 
tation 



T[V] = I W(r)^[T/(r)]exp{-7e[y(r)]}, 



(14) 



where, for a delta-correlated Gaussian potential with the 
correlator (||), 



(15) 



For large ui, the leading contribution to the sum ( p^ ) 
comes from transitions between the states Itpt) and 
with energies Et and i?b at the top and bottom of the 
Landau band, respectively. 



Et,h^ y"dry(r)|i^t,b(r)|' 



(16) 



To logarithmic accuracy, the conductivity is given by the 
solution of the variational problem of finding the optimal 
potential V(r) which minimizes the functional TZlV] and 
maximizes the matrix element of the transition subject 
to the constraint Et ~ E\j = fko, i.e., 

a{uj) oc max |exp [-7^[T^] -I- A {Et - E), - huj)] 

xKV-tlVl-lT/^b)!'}, (17) 

where A is a Lagrange multiplier. Variation with respect 
to V{r) gives the equation 

^ = A (l^tP - M + ^\nmVV\A)\' (18) 

(for brevity, we do not give the explicit form of the last 
term). 

We have analyzed the variational problem using a sim- 
ple and tractable direct variational method, and also by 
finding the maximum in Eq. ( p^ numerically. To see 
the qualitative features of the solution, we first discuss it 
ignoring the contribution of the matrix element. In this 
case the Lagrange multiplier A is given by the consistency 
equation. 



huj = Et-E^ = v^X J dr (IVtP - M) , (19) 
and then the conductivity ^% is 



I In 17(0') I = h Lu 



2v^ dr (l^tl'-MY 



(20) 



The solution ([T^) corresponds to a potential of the 
form of a well and a hump, far away from each other (cf. 
Fig. |l|). The potential is antisymmetric, the well and the 
hump have the same Gaussian shape [oc exp(— r^/2Z^), 
with r counted off from the corresponding extremum] 
and opposite signs. The wave functions V't and tpi, are 
localized at the hump and the well of V{r), respectively, 
and are given just by the most "localized" wave function 
of the lowest Landau level, namely that with zero angu- 
lar momentum, 7/'oo(r) oc exp(— r^/4Z^), centered at the 
appropriate potential extremum. The overlap of these 
wave functions is negligibly small, and Eqns. (|l9|), ( [20| ) 
give 



nLu = 2v'^XA, |lnCT(w)| 



with A = J c?r|-0t,b|^ — {4:TtP)~^. In this way we 
recover the expression (|l3|) for the conductivity tail. 
For higher Landau levels {N > 1), the wave functions 
have the formll3 t/^t^h oc exp{—iN(j)) cxp{—r^ 
in which case the corresponding constant An — 
(47rZ2)-i(27V)!/22^(Ar!)2. 

The prefactor in Eq. ([l^) prevents the well and the 
hump of from being too far away from each other. 
Nevertheless, the full variational equation ([ij) has a so- 
lution with an antisymmetric optimal potential Y(v) — 
—V[—v) and symmetric wave functions V't(r) — V'fc(— r); 
respectively, Et — —E\j = Tilo /2. To estimate the role of 
the overlap integral we used the direct variational method 
in which we sought the potential in the form V{y) — 
F(|r- rol) - F(|r-Hro|) with V{r) = exp(-rV2P). 
The distance 2ro separating the hump and the well was 
used as a variational parameter. Given the potential, 
one has to solve the Schrodinger equation, looking for 
the wave functions projected on the lowest Landau level. 
We took the functions 'i/'t.b in the simplest form of orthog- 
onal combinations of the zero-momentum wave functions 
centered close to ±ro (the positions were found using a 
variational procedure). The distance scales with fre- 
quency logarithmically. The overall asymptotic expres- 
sion for the exponent in a was the same as in Eq. (p^); 
the overlap integral gave only a prefactor, 

|(^t|VF|Vb)|' (?i7VM'ln(^/7)- 

[An extra w-dependent contribution to the overall pref- 
actor in a comes from the prefactor in the path inte- 
gral (p^). It actually increases with the increasing lj. 
However, an evaluation of this prefactor goes beyond the 
scope of this paper, and in some sense is superseded by 
the results obtained below with the method of moments.] 
To further check the accttracy of the asymptotic behav- 
ior of o'(aj), we maximizedllJ the functional in Eq. ( p^ nu- 
merically. We used the variational equation ([l^ ) to repre- 
sent the optimal potential as a bilinear combination of the 
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LLL wave functions V'Om(r) cx r™ ex-p{im(j)) exp(— 
with different magnetic quantum numbers m > 0, 



V(r) 



The corresponding eigenfunctions 'ipt,b were written as 
Unear combinations of the same functions ■0om(r). 

Both the exponent and the prefactor of the variational 
functional (|l7) calculated numerically become close to 
the result of the direct variational method for aj/7 3. 
The shape of the optimal potential found numerically for 
two values of u/j is illustrated in Fig. ^. 



a) 4 
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FIG. 3. Density plot of the optimal potential for uj — S-y 
(a) and tj = 87 (b). The distances are measured in units of 
the magnetic length /. 



B. Conductivity at small frequencies. 

An entirely different set of states defines the conductiv- 
ity for very small frequencies, <C 7. In this regime the 
constraint imposed by conservation of energy is not very 
restrictive, and it is the matrix element that determines 
relative contributions of different pairs of states. Close 
to the static limit the contribution to the conductivity 
increases dramatically with the typical size of a wave- 
function. As a result, a [Eq. (^)] is primarily determined 
by a narrow energy interval at the center of the Landau 
band where the states are nearly delocalized. The energy 
of the band center {E = 0) is a critical energy, similar to 



the critical value of the control parameter in the theory of 
classical percolation transition. At small deviations from 
the critical energy (parametrized by dimensionless energy 
£ = E/h"j) the correlation length diverges, ~ LjeJ|Zi|', 
where v = 2.33 ± 0.03 is the localization exponentfi'EHlS. 

Were it not for localization, propagation of a 
wavepacket in a random potential would be diffusive. 



(Ai?2(t)) ^ Dt. 



(21) 



Localization corrections are least important near the crit- 
ical energy, but even there theji modify the form of a 
wavepacket at very large timesllj. However, for not too 
large times the r.m.s. displacement, which is primarily 
determined by small momenta, retains the diffusive form. 
This can be used to find the conductivity at small fre- 
quencies. 

To this end, let us rewrite Eq. (foh 



2u^ 



lim SRe 

7 .5^+0 



dt e 



iujt~5t 



TroAi?2(t), (22) 



in terms of the squared displacement /S.R^{t) = 
[R(t) - R(0)]^ where R(t) = e'^^He^'"^ is the Heisen- 
berg operator of the guiding center. For an eigenstate 
\n) of the Hamiltonian (^) randomly chosen not too far 
from the critical energy, {n\/^B? {t)\n) has the diffusive 
form ( ^l|) at small enough but it eventually saturates 
at the distance of the order of the localization length . 
Replacing the trace by the integral over energy weighted 
with the (non-critical) density of states, we obtain the 
overall long-time (7t ^ 1) r.m.s. displacement 

TroAi?2(i) j de pih-fe) mm{Dt, $f ) 

cxZ2(^t)i-i/(2.)_ (23) 

This average is determined by the states with energies 
kl ^ (7^)^^^^"^; the integral rapidly converges outside 
this region. 

With asymptote (p3|), time integration in Eq. (|22 
gives: 



(24) 



Tlie.pSame result can be obtained from the scaling 
forniaa of the zero-temperature conductivity of the 
non-interacting system at a given chemical potential, 
which can be written as 



(25) 



where the dimensionless scaling function Qq{X) rapidly 
vanishes for X — s- and approaches a constant value for 
large X. Indeed, the conductivity for f3uj <C 1 can be 
written as a convolution of the scaling function (^) with 
the derivative of the Fermi distribution function 
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de 



(26) 



[cf. Eq. (||)]. For ksT ^ h^, all energies within the stripe 
of width Se ^ (lo/^)^ contribute equally, and in the limit 
w ^ we obtain Eq. (124) with the coefficient 



C = 16nfi 



dX 



GoiX). 



(27) 



Here we have assumed that e'^^e 
GoiX) = \im,^oGo{Xe^''ee/l^) 
verges both at zero and infinity. 



^ const for e ^ 0, and 
The integration con- 



C. Spectral moments 

Since the single-particle conductivity goes to zero both 
for cj ^ and for w ^ 7, its frequency dependence dis- 
plays a peak, with a maximum at a nonzero frequency 
w ~ 7. Such a peak is of central interest from the point 
of view of experiment, it does not occur in the SCBA. 
This peak was previously found and briefly discussed 
in our communication^. Here we present the results 
and provide some details of the full calculation of the 
low-frequency conductivity based on the method of spec- 
tral moments (MOM) . The advantageous feature of this 
method is that, instead of solving the full time-dependent 
problem of the electron motion in a random field (^ , one 
has to evaluate only equal-time correlation functions. 

We first calculate the spectral moments of the reduced 
conductivity (t{lo) [Eq. (M]. They are defined as 



1 

27r7 



duj {uj 



(28) 



For w, 7 ^ kBT/h, the states within the broadened Lan- 
dau level are equally populated, and the conductivity is 
an even function of frequency, a'(w) = a^—u). Therefore 
all odd moments vanish, M2k+i = 0. For even moments, 
we use the Hamiltonian (||) to obtain from Eqns. (11), 



(29) 



giq2fc+iR-l giq2ic+2R. 



The summation is performed over all qi , . . . , q2fe+2 • The 
commutators ( p9| ) can be evaluated recursively using 

[e*i^, e*i'^] = 2z sm(^h\ A q') e'^^+i') ^. (30) 

For Gaussian random potential, the disorder average 
in Eq. (29) can be computed by Wick's theorem. From 



Eqns. dFg) 



{VM = (t^IV^S) exp(-/2<7V2) Vq'.o, 



where S is the area of the system. Then, 

2 X k+2 



/ ;2 \ r 



C({q}) 



(I 



l2k+2 







■ qi A q2) sin ^— (qi + qa) A qs^ . . 
X sin (qi + q2 + • • • + q2fe) A q2fc+i j (31) 

where the sum is taken over all (2fc + 1)!! ways to choose 
pairs out of the set of 2fc + 2 variables, and 

C({q}) = '5(q.,i + qjj . . . 5{q,^^^ + q^.+J 
is the corresponding contraction function. 



a 



^ — )^ 



-X- 



^1 ^2 ^3 



% % % 



=X K= 



=X- 



-X- 



-X- 



\ 

-X- 



^2k-2 ^2k+2 

FIG. 4. Examples of diagrams. Dashed lines indicate which 
variables are paired, double lines represent an arbitrary inter- 
nal structure, (a) A symmetric diagram. Diagrams whose 
contribution is equal to zero: (b) disconnected diagrams, and 
(c) diagrams which vanish because the associated function is 
odd in q2fe+2. 



To classify different terms in the sum (|31|), it is conve- 
nient to depict the contraction procedure graphically as 
illustrated in Fig. |[ (These diagrams merely represent 
the contraction function C and should not be confused 
with Feynman diagrams for the Green's functions.) First, 
we note that the sequence of q^ paired in a diagram may 
be reversed, q^ — > q2/c+3-i, without changing the overall 
value of the diagram. The diagrams obtained by such 
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a reversal are equivalent, which reduces the computa- 
tion cost by a factor of two [this reduction does not oc- 
cur, however, for symmetric diagrams shown in Fig. 
Additional simplification comes from the fact that dis- 
connected diagrams (Fig. |b) and the dia grams with the 
structure shown in Fig. |jc are equal to zero. The number 
of diagrams of different sorts is given in Ref. |l^; for ex- 
ample, for the 14th moment there are 2027035 diagrams, 
out of which 5937 are symmetric, 318631 are discon- 
nected; the contribution of the diagrams with qi = — q2 
is « 72.559, whereas the contribution of all other dia- 
grams is ~ —8.809. 

Despite the reductions, the number of terms to be cal- 
culated remains very large for large k. Moreover, each 



term in Eq. ( ^T| ) is a sum of 2^'^ Gaussian integrals. Each 
integral can be calculated algebraically but at a high 
computational cost. To accelerate the calculation, we 
have devised an efficient numeric classification scheme, 
which sorts diagrams inexpensively into bins according 
to their approximate values calculated with double preci- 
sion. A representative diagram is evaluated algebraically 
for each bin. Finally, the diagrams are summed up with 
proper multiplicity, giving exact numerical values of the 
moments. The procedure is outlined in Appendix Cal- 
culating algebraically only non-equivalent Gaussian in- 
tegrals reduces computational time tremendously. This 
allowed us to evaluate the moments up to M14. For 
k = 0,1, ... ,7 we obtain 



M _ 1. 3. 443 25003 13608949709 298681273551508807 566602308094143977186611746328323669809 
^'^ ~ '8' 1152' 38400' 8941363200' 66698308912435200' 36033364452669289726755567308636160000' 
2589008911677049308284617052653287524724669331093372792412270459939701 

40611974008223423608381355617240666314144290787406293503186042880000 ^ -* 



and the corresponding approximate values, 
M2k ~ 1; 0.375; 0.385; 0.651; 1.522; 4.478; 15.72; 63.75 . 



D. Reconstruction of frequency-dependence 



Since the conductivity is asymptotically Gaussian, one 
is tempted to restore a{Lu) from the moments M„ in a 
standard way, writing an expansion in Hermite polyno- 
mials cr(7x) = Bn Hn{\/2x) exp{—2x^) . The coeffi- 
cients Bn can be expressed recursively in terms of the 
moments M^, k < n. However, for the moments (|32]), 
this expansion does not converge rapidly, see Fig. |[ This 
is consistent with nonanalyticity of the conductivity at 
w = 0. 

Given the exponent ^ in Eq. (|2j), a much more rapidly 
convergent expansion can be constructed in terms of a 
different set of orthogonal polynomials. Specifically, with 
Eqns. (|l^), (|2^), we write the conductivity at the lowest 
Landau level as 

a{uj) = x^'G{x) exp(-2x2), x = \uj\/-f. (33) 

The function G{x) (x > 0) can be expanded in Laguerre 

polynomials L^l^ ^•'^^(2a;^), which are orthogonal for the 
weighting factor in Eq. ( p3| ) . It is important that the ex- 
pansion coefficients can also be recursively restored from 
the moments M2k, k < n. 
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FIG. 5. Approximating d-{u>) with Hermite polynomials. 
With the increasing number of moments M2k, a minimum 
of the conductivity is consistently developing at = 0, and 
the expansion does not show fast convergence for small uj/"/. 

For the presently accepted value of the localization ex- 
ponent v ~ 2.33, the value of the conductivity exponent 
fj, is 0.215. The expansion for G converges rapidly for /x 
between 0.19 and 0.28, whereas outside this region the 
convergence deteriorates, as illustrated in Fig. ||. This 
could be considered as an indirect indication of the con- 
sistency of our approach. 

The resulting conductivity calculated with /i — 0.215 
is shown in Fig. |^. The estimated deviation from the 
obtained curve due to the finite number of moments and 
also to the uncertainty in the value of fj, (its effect is dis- 
cussed in Ref. is smaller then the width of the line. 
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FIG. 6. The expansion of the prefactor G (|3^) in Laguerre 
polynomials ^^^'^{2x^) for different values of the exponent 
/i, depending on the total number of moments M2k- The ex- 
pansion converges rapidly for ^ between 0.19 and 0.28 and 
much slower outside this interval. 



III. SINGLE-ELECTRON CYCLOTRON 
RESONANCE 



Resonant single-electron absorption near cyclotron fre- 
quency {uj « LOc) is determined by the correlation func- 
tion of the Landau level raising and lowering operators 
p± defined in Appendix 



ne2 r 
2m J_, 



dt e*'^* 



(p_(i)p+(0)), 



(34) 



where we used Eq. ^ and assumed that 
eiijp{hujc/kBT) ^ 1, in which case only the lowest Lan- 
dau level is occupied for small densities. If the disorder 
is weak, 7 <C Wc, it only weakly mixes different Landau 
levels, the primary effect being lifting the degeneracy of 
each level. Then, the problem of optically-induced tran- 
sitions between different Landau levels resembles that of 
transitions between degenerate electronic terms of im- 
purities in solids in the presence of the electrpH-phonon 
coupling which gives rise to Jahn- Teller effectEHl. One of 
the effective methods of the theory of absorption spectra 
of Jahn- Teller centers is the method of moments. 

The MOM formalism can be transferred to the case of 
inter-Landau-level transitions. For ^ fc^T, quenched 
disorder can be described in the same way as scattering 
by thermally excited phonons. The major difference is 
infinite level degeneracy. 

In the neglect of disorder-induced scattering between 
Landau levels, one should keep only the part Hd of the 
disorder potential V{r), which is diagonal in the Landau 
level representation, 



N q Af 



(35) 



where Vq is defined by Eq. (g) and Pjy — P% is the opera- 
tor of projection to the A^th Landau level, as in Eq. (|A2|). 
With the Hamiltonian (|35|), oscillations at the cyclotron 
frequency can be singled out in Eq. (p3). 



Then, from Eq. (jsj), we can write 

„2 



-iHdt/h 



O-c(w) 



ne 
2m7 



(36) 



(37) 



where 
o'c(w) = 7 



Tro [e^Hdt/hp_^-^Hdt/hp^-j (38) 



is the reduced conductivity, and Alo = — aj^ is the 
frequency detuning, |Aa;| <C uJc- 

The major difference of Eq. ( ^8| ) from its counterpart 
( p!o| ) for the low-frequency conductivity is that the Hamil- 
tonians for direct and inverse time propagation (corre- 
sponding to the factors e**^"*') are now different, which 
is again familiar from the theory of impurity absorption 
spectra. The reduced cyclotron conductivity can be con- 
veniently written in a form conventional for this theory 
by introducing the "perturbation" Hamiltonian 



(39) 



In the interaction representation, CTc can be then simply 
expressed in terms of a time-ordered exponential. 



aduj) = 7 



dt e 



iAu t 



Tro 



T,exp(- dT(5i/d(r))J. 

(40) 

Here, time dependence of the operator 5Hd, 



(41) 



is generated by the disorder Hamiltonian projected on 
the LLL, which is given by Eq. (|^) of the previous sec- 
tion. 

We can now define the spectral moments of the cy- 
clotron peak as 



1 



duj 



27^7 7-c 
Using Eq. (EO) we write 



(42) 



8 



t=o 



(43) 



We note that, similar to the case of the peak of low- 
frequency conductivity discussed in the previous section, 
we are calculating here the moments of the cyclotron 
peak only, whereas the small (cx ^ /oJc) background from 
the correlators neglected in obtaining Eq. ( p^ ) is pro- 
jected away, as are also the peaks of Uxx{'^) at a; ss nojc 
with n 7^ 1. 



A. Tails of the cyclotron resonance peak 

As in the previous section, let us first discuss the 
asymptotic form of the cyclotron peak comparatively far 
from resonance, |Ati;| ^ 7 (yet jAcjl ^ lOc)- If we in- 
troduce the exact eigenstates of the Hamiltonian (|35| ) 
for the lowest |0, m) and the first excited |1,to) Landau 
levels, with energies and respectively, the ex- 

pression (|8|) for the reduced conductivity can be written 
in the form 



with different Hamiltonians for Et and Eb\. If we ne- 
glect the dependence of the transition matrix element on 
V{y), then for the optimal potential we obtain an equa- 
tion similar to Eq. (45). However, in contrast to Ref. 
we permit the centers of the wave functions $0 and $1 
to be shifted with respect to each other. 

A remarkable feature of this simplified variational 
problem is that, in the neglect of overlapping of the dis- 
placed wave functions, the same value of the variational 
functional [except for the overlap term] is obtained for 
the trial wave functions of the first Landau level with the 
magnetic quantum numbers -1 or 0, i.e., '01,-1 or V'l.Oi 
or for any their linear combination. 

For a displacement R between the centers of the hump 
and well of the optimal potential, the transition matrix 
element is KV'ilp+l'^o)! ~ exp(— i?^/4/^). The optimal 



distance ~ ln[{uj 



' /"f ] is found by maximiz- 



ing the expression with the matrix element present. As in 
the case of the low-frequency conductivity, this distance 
increases as the frequency is tuned away from resonance. 
The variational result for the conductivity tail is 



CTc(cj) oc exp ^ - -^{'^ - ^c)'^Y 



(46) 



This tail is much broader, with the exponent reduced by 
a factor of 3, compared to the result of Ref. n4. 



(44) 

As for the low-frequency conductivity considered in the 
previous section, the conductivity tail is determined by 
large optimal fluctuations of the disorder potential. 

The problem of the optimal potential for cyclotron res.- 
onance was previously considered by loffe and LarkinliJ. 
They used an ansatz of a rotationally-symmetric optimal 
potential 



277^1 I $1 



(45) 



where $0 — "00,0 (r), 'I'l — "01,-1 (r) are the functions 
of the lowest and first excited Landau levels centered 
at the same origin, with magnetic quantum numbers 
and —1, respectively. This resulted in the asymp- 
totic form of the cyclotron resonance absorption peak 
(Tc oc exp(— 8 Aa;^/7^), for the range <C ksT. 

We argue that the transition probability between the 



states with energy separation Em^ — E, 



AO) _ 



Aw 



-I- LUc is 



exponentially increased if the cyclotron orbit centers of 
these states are permitted to shift with respect to each 
other. This happens despite the associated decrease of 
the overlap integral. 

The calculation of the tails of the cyclotron reson ance 
absorption peak is very similar to that in Sec. [I A . We 
begin by writing the averaging in terms of a functional 
integral (^^ , with the energy conservation taken into ac- 
count using a Lagrange multiplier [as in Eq. (17) but 



B. The center of the cyclotron absorption peak 

Generally, we do not expect CTc((j-') to display a non- 
analytic dip at the center of the cyclotron absorption 
peak. Indeed, the power-law singularity (|2^ ) of the low- 
frequency conductivity can be associated with quantum 
interference which leads to eventual localization of all 
states except for one (or maybe a few) at the band cen- 
ter. The expression (^) for the cyclotron resonance ab- 
sorption has a structure which differs from that for the 
low-frequency conductivity. In particular, it contains an 
extra phase factor from the Hamiltonian SHd ( p9| ) which 
represents the difference in the random potential expe- 
rienced by an electron at the two Landau levels. This 
phase factor should give rise to an exponential damping, 
and related suppression of the interference effects at long 
times. Consequently, the conductivity is expected to be 
smooth near lUc- 

Another way to see this is based on the following ar- 
guments. The suppression of the low-frequency conduc- 
tivity for w — > may be attributed to level repulsion 
between overlapping localized states. This repulsion is 
comparatively small for states of relatively large radii, 
with energies close to the band center. Indeed, only such 
states contr ibute to the low-frequency conductivity, as we 
saw in Sec. 11 B. On the other hand, resonant cyclotron 



absorption is due to transitions between different Lan- 
dau levels. Although the central part of the absorption 
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peak is formed by transitions between strongly overlap- 
ping states, the involved states are eigenstates of different 
Hamiltonians, with random parts and h'^^\ Their 
wave functions have different spatial structures and their 
energies are essentially uncorrelated, except for states 
deep in the tails of the Landau levels. Consequently, we 
expect no suppression of transitions at frequencies close 
to ujc- This argument is in agreement with the results of 
the method of moments presented in the next subsection 
and Fig. I 



C. Spectral moments of the cyclotron peak 

We will now calculate the spectral moments (^). Be- 
cause all states of the lowest Landau level are equally 
populated, the reduced conductivity ( |40| ) is symmetric 
with respect to Wc, (^d'^c + Aoj) = crd^c — Acj). 
Therefore all odd moments vanish, M^j^j^^ = 0. The 
structure of the expression for even moments whi ch fol- 
lows from Eq. (^) is similar to that of Eq. (Bl). The 
main difference is that the prefactor now is a complicated 
polynomial, a combination of products of terms which are 
linear in the squared wave numbers q\. 

This strongly complicated the numerical procedure. In 
particular, we failed to find any symmetries to reduce 
the computational overhead, and graphical representa- 
tions were of little help. We were also unable to catego- 
rize different terms as described in Appendix ^ for the 
low- frequency conductivity. InsteacLwe developed the 
computer algebra package Gauss Intc3 for Mathematica, 
capable of handling the integration of high dimensional 
Gaussian integrals in a manageable time frame, and used 
the brute-force approach calculating all terms analyti- 
cally. For fc = 0, 1, . . . , 5 we obtained: 



1 37 52043 4750893001499 



(47) 



2' 64' 55296' 2488320000000' 
29694054188353275207831950716496054687 

6480696333914117611721116876800000000 ' 
and the corresponding approximate values 

M|fc « 1.000; 0.500; 0.578; 0.941; 1.909; 4.582. 



The values of the moments with k — 0,1,2 were also 
independently confirmed analytically. 



The coefficients Bn were recursively expressed from the 
calculated moments (0). We discovered that although 
the convergence is fast far from the center of the peak, 
it is noticeably slower close to the center [we emphasize, 
however, that we reached convergence, in contrast to the 
similar expansion for the low-frequency conductivity in 
Fig. H, where the convergence was not reached for 14 mo- 
ments] . Within this approach, the number of calculated 
moments is apparently insufficient for restoring the entire 
function dc with desired accuracy. The corresponding re- 
sult is shown with dashed line in Fig. |^). 

Much faster convergence was achieved when ^d'-^) was 
restored using a continued fraction expansion. We ap- 
plied an algorithm similar to that used to reconstruct the 
LLL density of states from its frequency, moments for an 
arbitrary correlated random potentialEj. The steps in- 
volved in this process are summarized in Appendix |^. 
As one can see from Fig. |^, the convergence is very fast. 
The resulting shape of the cyclotron absorption peak is 
shown in Fig. ||. We believe that the deviation from the 
exact value is within the width of the curve. 




-1 



(co-coj/y 
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D. Reconstruction of frequency-dependence 

As a first step, we reconstructed (Tc(t^) using a stan- 
dard expansion in Hermite polynomials, 

adix) = ^ S„ i/„(v/873x) exp(-8a;V3). 



FIG. 7. Approximating Oc with continued fractions, (a) 
The convergence with increasing number of moments Mjj. 
is extremely fast: the curves lie on top of each other al- 
ready for k — 3,4, 5. (b) A comparison between the con- 
tinued-fraction (solid line) and the Hermite polynomial ap- 
proximation (dashed line) for fc = 5. Convergence is much 
faster with continued fractions. 
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IV. CONCLUSION 

In conclusion, we have analyzed the low-frequency 
single-electron magnetoconductivity and cyclotron reso- 
nance absorption of a nondegenerate 2D electron system 
in a quantizing magnetic field. We considered the exper- 
imentally important parameter range where the width of 
the Landau levels is less than temperature, so that all 
states within the lowest Landau level are equally popu- 
lated. In this range, by combining the ideas of the scaling 
theory of the IQHE, the method of optimal fluctuation, 
and the method of spectral moments, we obtained highly 
accurate numerical results throughout the frequency do- 
main where the conductivity displays peaks. 

We found that, in contrast to tiie prediction of the 
SCBA or other mean-field theoriesE^, the low-frequency 
conductivity displays a peak at a nonzero frequency, as 
shown in Fig. |^. For short-range disorder, the position 
of the peak is given by 



0.267. 



(48) 



For a; — > 0, the single-electron conductivity displays a 
universal power-law dispersion Uxx oc lo^ , which is re- 
lated to the scaling behavior of the localization length 
as a function of the distance in energy from the center 
of the disorder-broadened Landau level. On the other 
hand, the peak of the cyclotron resonance does not dis- 
play such singular behavior and is not shifted away from 
Wc, as seen from Fig. ^ Both peaks have Gaussian tails, 
with different exponents [see Eqns. (|l^), (^)]. 

Experimentally, it is more feasible to investigate the 
magnetoconductivity at a given nonzero frequency lo as 
a function of the external magnetic field B. The corre- 
sponding representation of our results is given in Fig. ^ 
for the scaled conductivity (T*(i3;w), 



(7* = (T*(B; u) 



B 



1/2 ~ 



ct(7)^ 



taken into account, this dependence becomes steeper, 
with B^f^a^ cx B^''/^ as illustrated in Fig. |[ 

Within the single-electron approximation, the restric- 
tion on the magnetic field from above is imposed by 
the condition Ti^ <C fcsT, which is equivalent to tOc ^ 
To (kBT/h)'^, for the short-range disorder potential. This 
inequality can be fulfilled simultaneously with B ^ 
Bq{lu) provided huj <C fcsT. The restriction on the mag- 
netic field from below necessary for the system to be 
in the lowest Landau level, hwc 3> ksT, can hold for 
B ^ Bq{uj) and fiuj <^ ksT provided ksT ^ ^Tq^. 

The latter inequality is often fulfilled for electrons on 
helium surface. In this case, for T < 0.9 K the random 
potential V(r) is due mostly to capillary waves, ripplons. 
It has a small correlation length and is quasistatic. For 
electron densities n w 0.5 x 10* cm~^ and T = 0.7 K, 
the value of tq is as big as w 2 x 10^* s (see Refs. p|p6|). 
For lower T, the mobility which corresponds to the ef- 
fective To « 10~^ s has been observed by Shirahama et 
al.cli. The condition fiTQ^ <C fc^T can thus be easily 
met. Therefore the results of the present paper fully ap- 
ply to electrons on helium as long as one can use the 
single-electron approximation. 



(49) 
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where the scaling factor (7(7) « 1.08, and the scaling field 
Bq(uj) is defined by the equation j{Bq) — uj. The mag- 
netoconductivity axx is related to (T*(i?;w) by a factor 
which is independent of B (but depends on w), 



(7(7) h ne^ 



47r ksT mojTo 



cr* {B;uj). 



(50) 



Here, t^^ = mv-^/h^ is the rate of electron scattering 
by the short-range potential (^) in the absence of the 
magnetic field. The frequency-dependent scaling field 
in Eq. (^9|) is related to uj and tq by the expression 
Bo{lu) — TrmcTQUj'^ /2e. 

In the self-consistent Born approximation, the func- 
tion fT*(i?;w) decays with the increasing magnetic field 
as B~^''^, for B > Bo{uj). With the localization effects 



FIG. 8. Reduced ac magnetoconductivity cr* ( p9| ) at a 
nonzero frequency uj as a function of tlie reduced magnetic 
field B/Bo{lu) oc B uj~^ . In order to demonstrate the anoma- 
lous single-electron behavior, a, is also plotted with an extra 
factor {B / BqY^^ . For large B, the single-electron conductiv- 
ity displays scaling behavior, B^^^a, oc B~^^^ . 

The most restrictive limitation on the present theory is 
imposed by many-electron effects. Even where the elec- 
trons do not form a Wigner crystal and the 2DES is non- 
degenerate, thesa ^effects can determine magnetotrans- 
port phenomenaQ'ES. However, for sufficiently strong 
magnetic fields, the effective coupling to a short-range 
disorder potential becomes in some sense stronger than 
the electron-electron interaction. Therefore the single- 
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electron approximation describes certain features of the 
magnetoconductivity for strong fields. In particular, al- 
though the fuILmany-electron static magnetoconductiv- 
ity is nonzerc£3, axxi^^) can still display a peak at the 
frequency (p8|). Observation of this peak and/or its coun- 
terpart in the magnetic field dependence of the weighted 
ac conductivity B^^'^axxito) (cf. Fig. ||) would be a clear 
demonstration of single-electron localization effects in 
quantizing magnetic fields. 
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valuable discussions. Work at MSU was supported by the 
NSF through grant no. PHY-00071059 and by the Center 
for Fundamental Materials Research. L.P. was supported 
in part by the DOE grant DE-FG02-90ER40542. 



APPENDIX A: GUIDING CENTER FORMALISM 

Here we remind the basic expressions of the guid- 
ing center formalism needed to derive simplified expres- 
sions for the low-frequency conductivity at a given Lan- 
dau level. Time evolution of the momentum operator 
p — —ih'V + (|e|/c)A in Eq. (^) is determined by the 
Hamiltonian Hq + V{r) where V{r) is the scattering po- 
tential and 



Ho = hujc (p+p- + ^ 



Here, p± are the Landau level raising and lowering oper- 
ators 



P± 



{Px T ipy)l\/2nmujc, [p-,p+\^l. 



The choice of signs corresponds to i?2 = — |B| < 0. 
The guiding center coordinates are defined as 



R = (X, r) : X 



Py 



Px 



They commute with the momentum operator, [R, p] = 0, 
but not with each other, [X, Y] = —iP. To express the 
disorder potential in terms of these variables we expand 
it in Fourier series, V{y) = Kjexp(iqr), which gives 



Q^i^/i Jq-p+ „-iq+p- 



(Al) 



where q± = {qx T 'iQy)/V2- Introducing the operators of 
projection to the nth Landau level Pn, we can write the 
projected disorder potential as 

V^^Hr) = Pn V{v)Pn = h-,PNj2 l^r^e^'^'', (A2) 

q 

V-W = (Vjh^) L^iq' I'm cxp(-g2;2/4), (A3) 
where Lfq{z) is the iVth Laguerre polynomial. 



Using the commutation relation [e*'i"',p±] = lq±e'^^'^, 
we can evaluate the Fourier-transformed current-current 
correlator in Eq. (0) as 



{px{t)px{0))^ 



where we use the notation 







1 - 


(5)1 







{px{t)pxm. 



dte'^\px{t)pxm- 



The correlators in the r.h.s. of Eq. (A4) can be iterated 
similarly, and the whole current-current correlator can be 
expressed in terms of the correlators of the derivatives of 
the potential V (cf. Ref. |l3|). 

In the adiabatic limit, uj <^ lOc, the first term in the 
r.h.s. of Eq. ( |A4| ) vanishes, and from this equation or di- 
rectly from the equation of motion, for times slow com- 
pared to wj^, we can identify 



1 d 



Pxit) j^Vm] = —J^qyV^e^^'^^'K 



The rightmost part of this expression is proportional to 
the time derivative of the guiding center coordinate, 

x(t) = ~{i:'in)Y^q,Vqe^^^'^ 



[in the right-hand sides of the last two equations one 
should keep only the terms which do not contain fast- 
oscillating factors exp(i7Va;ci) with iV 7^ 0]. We can then 
further identify the projected part of p with mR, which 
gives Eqns. (|lg) and (0). 



APPENDIX B: CALCULATION OF FREQUENCY 
MOMENTS 

Here we describe an efficient numerical procedure used 
to obtain exact expressions of high-order diagrams in the 
expansion (|3l|). 

We begin by rewriting Eq. ( |3l| ) as a sum of exponen- 
tials, 



M2fe = -2/2 



^ C({q}) {b} 



X y^qi . . . rfq2fe+2 C({q}) (qiq2/c+2) 



2 2fe+2 



(Bl) 
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The inner sum is taken over all binary sequences b = 
(61, 62, • ■ ■ , b2k), with bi = 0, 1. They label possible com- 
binations of signs which arise from writing the sines 
in Eq. ( ^T]) in terms of exponentials. The quantity 
cr({b}) = J^i^i- The antisymmetric (2fc + 2)-dimensional 
square matrix B^^^ has the following structure: 



/ 



^{b} 





-ci 

~C2 



-C2k 

\ 



Cl C2 
C2 
-C2 



-C2k 



-C2k 



C2k 
C2k 



C2k 









0/ 



(B2) 



where Ci = (—1)''*. 

Because of the (5-functions in the contraction function 
C({q}), integration in (Bl) has to be performed over fc + 1 



independent wave vectors. Up to a prefactor (qi q2fe+2), 
the integrand is an exponential of the quadratic form 
^ q^Ay q^ , where i,j = l,....fc+ 1. The ma- 
trix elements Aij are themselves 2x2 matrices, Aij — 



aij&y, where ay is the Pauli matrix, and aij 



0, ±1. 



For a given contraction C({q}) and a given vector b in 



(Bl), i.e., for the corresponding matrix , the Gaussian 



integrals can be evaluated exactly, giving 



k+l k+1 



1 + A2 



(B3) 



where iXm are the eigenvalues of the antisymmetric ma- 
trices a„m, and u„i,n SLie the components of the corre- 
sponding eigenvectors. The subscript x takes on the value 
a; = fc-|-l if qi and q2fc-i-2 are independent variables in the 
pairing procedure, whereas for qi = — q2fe+2 we should set 
X = 1 and additionally multiply Eq. (B3) by (—1). 



Given the large number of terms and the computa- 
tional price of calculating the exact values in each term, 
we did not calculate each integral exactly. Instead, the 
value of a given term (specified by the choice of contrac- 
tion C({q}) and the binary sequence b), was calculated 
numerically with double precision. With integer-valued 
matrices a^ , the number of different values was not ex- 
ceedingly large, and we used the obtained approximate 
values (B3) to assign each term to an equivalence class. 
The individual weights can be positive or negative, de- 
pending on the parity of (T({b}). This procedure was 
used instead of much more tedious manual classification 
of high-order diagrams. (In addition, we checked for sev- 
eral values of k that the weights corresponding to special 
diagrams in Fig. ^ are indeed equal to zero). 

After the classification of diagrams was completed, the 
exact value of the integral in each class was obtained by 
taking a representative contracted matrix aij and cal- 
culating the Gaussian integral I[a] algebraically. These 
integrals have rational values; the final answer for M2k 



was obtained as a weighted sum of these rational num- 
bers with their respective bin weights. As a test, we 
compared the results for /c = 0, 1, 2 with explicit analytic 
calculation. 



APPENDIX C: CONTINUED FRACTION 
EXPANSION 

The Stieltjes transform of the conductivity ac is de- 
fined by 



duj 



(7c(w -I- LOc) 



z — iuj/"f 

while the inverse transformation has the form 



3?e z > 0, (Cl) 



adio + ujc) = 2 lim ^c[R{e 



The function R is related to the moments 
expression 



i?(z) = ^z'^M,^ 



,-fe-i 



(C2) 
by the 

(C3) 



k=0 



We now construct an approximation for (Cl) which 
applies for an even function a{uj + l Uc) — a{—Lu + ujc), 
allows for the Gaussian asymptotics (nw 



lim —liiac{Lo + uJc)^~Tr- 



LU — >±oo CO 



2a 



16' 



and requires aabi-.a finite number of moments. 

It is knownnj'E2l that an odd function R{z) can be ex- 
panded into a Jacobi-type continued fraction. 



Riz) 



K 



A, > 0, 



where we use the notation 
'A. 



K — 



(C5) 



(C6) 



z + 



A, 



z+ ■■ 

for the continued fraction with coefficients Aj and 
variable z. The first J continued- fraction coefficients 
Ai, . . . , A J are obtained from the normalized moments 



M2, ■ ■ ■ , M^j by expan ding the power series (C3) into 



the contiuiied fraction (C5) using an efficient recursive 
algorithmcj. Having obtained only a finite number of 
coefficients A,-, we need to estimate the remainins-enes. 



Fortunately, the asymptotic behavior (C4) implieS the 
following asymptotically linear growth for the continued- 
fraction coefficients 



13 



lim — - — a. 

j^ca j 



(C7) 



Therefore, if the first J coefficients Ai, . . . , Aj are found, 
one can then construct an approximation R^'^\z) to R{z) 
by finearly continuing Aj for j > J, 



R^'\z) 



K 



A 



(./) 



where 



for j < J 
A,7 + - J) for i > J 



(C8) 



(C9) 



A continuous fraction with a finearly increasing coeffi- 
cient can be written in terms of the Whittaker parabofic 
cyfinder function D^, 



i/2D_^/„(a-i/2^)' 



which is vafid if a > , /? + a > and 5fte 2 > 0, so that 
we can write Eq. (C8) as 



1 



z + 



Ai 



(CIO) 



z + 



z + AjT{Aj,a,z) 



Applying the inversion formula ( |C2| ) immediately gives 
the restored cyclotron resonance absorption (Tc(w) as 
shown in Fig. M. 
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